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LIFTING SURFACE THEORY FOR A HELICOPTER ROTOR IN FORWARO FLIGHT 
Harry L. Runyan and Hsiang Tai 
Abstract 

A lifting surface theory has been developed for a helicopter rotor in 
forward flight for incompressib 1 e flow. The method utilized the concept of the 
linearized acceleration potential and makes use of the vortex lattice 
procedures. Results in terms of lift coefficient slope for several forward 
flight conditions are given. 

1. INTRODUCTION 

Rotating lifting surfaces are an integral part of the propulsive unit of 
every aeronautical and nautical vehicle, from the compressor and turbine bladc-s 
of jet engines, the pumps for rocket engines, to propeller and helicopter 
rotors. The aerodynamics of these rotating elements has been under extensive 
study since the advent of the airplane and with a combination of experimental 
and analytical approaches, succcessful designs have been achieved. In many 
cases, two-dimensional theory has been used, usually modified by an assumed 
spanwise distribution, and inflow angles. A practical lifting surface theory 
has not been developed for the compressible rotating case for the general case 
of a helicopter in forward flights. It is the purpose of this report to present 
a lifting surface method within the limits of linearized theory. 

The method is based on the concept of the acceleration potential, 
originally introduced by Hans Kussner (ref. 1). The method was first applied to 
an oscillating wing in uniform translatory motion by one of the authors 
including compressible flow in 1957 (ref. 2). The acceleration potential 
approach has now become standard in the determination of the unsteady 
aerodynamics forces for flutter studies. 
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The first use of the acceleration potential approacr. for a rotating system 
was made in a paper hy Hanaok3 (ref. 3) for the Incompressible case In a paper 
directed at application to the marine propeller. The acceleration potential has 
been used in the past in studying the oropeller noise problem, but in all of 
these noise propagation cases the problem was specialized early In the 
analytical development to the so-called far-field case usually with the 
stationary observer, whereas the lifting surface theory is essentially concerned 
with the details of the near-field case for a co-moving observer as well as the 
satisfaction of certain boundary conditions. Hat, (ref. 4), has derived a 
general expression for an acceleration doublet for any motion. His development 
has been used as a basis for this investigation. The procedure developed here 
involves the precise numerical integration over the surface of the rotor. It is 
hoped that the method lays fundamental ground work for propeller and helicopter 
unsteady, three dimensional compressible aerodynamic theory. 

With regard to the contents of the report, the next section contains a 
brief derivation of the fundamental equations, including a discussion of some 
implications of the equations. The third section contains a description of the 
method of solution. Finally, the results of some preliminary calculation for 
the incompressible case are given. 

?. BASIC FPIIATIONS 

The fundamental equation required for the solution of the lifting surface 
problem relates the downwash velocity (velocity normal to the lifting surface) 
to the unknown pressures on the rotor blade. 
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The acceleration potential method will be used In the analysis. The 
acceleration potential, ♦, is defined as 


P - p. = -Pt (2.1) 

Rernoulli's equation for pressure in terms of the velocity potential, $, is 
defined as 


P - P £ 




( 2 . 2 ) 


where the operation represents the substantial derivative. For a wing 


translating in the x-direction with velocity II, the equation can be written 

p - p. - * "£> 

or 

"ft < 2 - 3 > 

For wings oscillating in harmonic motion such as in flutter, an harmonic 
substitution can be made to obtain 


^ = -( i to 



Solving for $ gives 


-itox 
e T" 


x 

/ 


i uiX 

*(i,y,z)e U dx 


(2.4) 


Since the motion of a rotating helicopter blade is non-harmonic, this 
procedure is not practical, thus the relationship is 


* 


P<f> 
= W 


Equation (2.5) can be integrated with respect to t. to give 

= / t ipdt 


(2.F) 


(2.fi) 
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’inally, the downwash is derived 

w ■ In ■ -ar ! * dt (2 - 7) 

-00 

3. ARBITRARY MOTION Of AN OSCILLATING DO H RLET IN COMPRESSIBLE FLOW 

The expression for an oscillating doublet moving In an arbitrary manner In 
compressible flow was derived by Oat (ref. -5). A concise derivation of the 
doublet expression is given below. 

The expression for an accelerating source is 



4„|? - P o (t)| 


1 


q ' i V T > 

J 0 (tH* - ? 0 (t)] 

a l p : 


(3.1) 


where £ 0 (t) the vector position of the source at time t, P is the position 
of the field point, t/ 0 ( t ) is the velocity of the source point at time t, 

"a" is the speed of sound and q is the strength of the source. 

An auxiliary equation which relates the time interval t - t between the 

time the source emits a disturbance at t and the time at which it is received 

at the field point at t, to the distance between the two points is given below 


\ [ - Vo 

a 


(3.2) 


A doublet may he formed by placing a source and a sink in close proximity 
and allowing the distance to approach zero, while at the same time maintaining 
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the product of the distance and the source strength finite. This result, may 
also be accomplished by performing a differentiation in the proper direction. 


Hat applies this second technique by first defining a vector distance n 



(3.3) 


where £ is a parameter and n 0 is the unit vector which is in the desired 


direction. Substituting D for P - P 0 in equation (3.1) gives 



(3.4) 


The doublet is now obtained by applying the derivative with respect to £ and 


letting £ ♦ 0 or 




► 


ORIGINAL PAGE IS 
OF POOR QUALITY 


6 



(3.6) 


This expression checks an expression derived in a different manner in ref. 7. 

4. INCOMPRESSIBLE PLOW. - The incompressible case may be obtained by 
letting a °° to get 


* 


D 


• % 
4ir | D | 3 


(4.1) 


the downwash equation may now be written 


w 



— oo 


q( r n.t) 



3(n • D)(n 0 • D) 

iSi* 


d T 


(4.2) 


The integral equation relating the known downwash, w(P,t) to the unknown 
pressures, q(r Q ,T), is given by the following equation for incompressible 
flow 


w (p,t) = 


1 

4 Tip Chord 


/ 


Span 


/ f l q(r 


o’ 


t) 




dc dr dt (4.3) 

-> c 0 

|D| 5 

In eq. (4.3), D = P(r,t) - p 0 ( r 0 »T) denotes the distance between the position 
vectors of the downwash and doublet points relative to a fixed coordinate 
system, where r, r Q are measured radially along the span; t,x are the 
independent time variables of the dowrwash and doublet point respectively. 
q(r , T ) is the lift strength per unit length of span; p is the density of 
the fluid, n and n Q are the unit normal vectors perpendicular to the 
surface at the downwash and doublet point, respectively. 
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The final downwash equation Involves a double integration over the surface 
of the blade which has proved to be a stumbling block to the efficient 
application of lifting surface theories. Most of the research directed towards 
the application of the kernel function for the planar wing In uniform motion has 
been aimed at devising feasible and adequate means for accomplishing this 
integration. Of the several techniques avallahle, the one selected for 
application to the rotor case is termed the doublet-lattice method, a principal 
assumption in the method is the placement of line loads at the 1/A chord of the 
wing and the satisfaction of *he downwash condition at the 3/4 chord location of 
each element. This would make eq. (4.3) tractable, namely by yielding a set of 
linear equations for q^ . The total lift of the blade at the time t, or 
equivalently at certain azimuthal position will be £qi(t). 
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Integration over t, it is tantamount to considering *.11 positions successively 
occupied by each lifting element in the past, which amounts to scanning the 
whole wake. The blade has the cord c end length Rt - Rs» ^s being the 
root of the blade, Rt is the maximum length of the blade. Suppose at the 
time origin, the blade momentarily coincide with the coordinate system along the 
x-axis as shown in the sketch. 
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The blade executes a counterclockwise rotation with angular velocity 
ft while novlng with velocity U along the negative x direction and W along 
the positive z direction. Since the quarter line method has been adopted, the 
doublet point lies 1/4 c ahead and downwash point lies 1/4 c aft of the 
nidchord. The positions of the doublet point as well as the downwash point can 
be established as follows. The Cartesian components of the doublet position are 


x 0 = -Ut + r 0 cos(ftx) - c/4 sin(ftx)cosa 0 
y 0 = r 0 sin(ftx)+ c/4 cos(nx)cosa 0 

z 0 = Wx + (c/4) sin a 0 (4.4) 

With the substitution of c ♦ -c, r 0 ♦ r, x ♦ t the position of the downwash 
point is given by 

x = -Ut + r cos(ftt)+ c/4 sin(ftt)cos a 

y = r sin(ftt)- c/4 cos(ftt) cos « (4.5) 

z = Wt - c/4 sin a 

In eqs (4.4) and (4.5), the angle a,a 0 are the twist angle of the blade 


defined below 
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tan ° " U sin(ftt)+ rn 


10 


w 

ta " °o ' II slnlnt)* r o i) 


(«.«) 


This means that the twist argle is determined by the velocity vector at the time 

t (i.e., the final time) which implies that the angle of twist is constant in 

time, and is only a function c* radial distance. Therefore, the reference plane 
(defined hy the doublets and downwash points) is a twisted surface. Only when 
the doublet (hereafter denoted as d.h.) and the downwash point (hereafter 
denoted as d.w.) possess the same radial distance, r = r 0 , then the d.b. and 

d.w. lie jn a plane (i.e., n jj n 0 at t = t). 

From eq. (4.4) the doublet velocity can be computed, namely the time 

derivative of the position vectors. 


♦ • ♦ 


V = xi + yj + zk 
o o ■' o° o 


(*. 7 ) 


where 


x o = " r o n s ^ n (^ T ) _ c /4fl cos(fli)cos a Q 


y o = r Q 0 cos Ot - cft/4 sin(Oi)cos a Q 


To determine n Q as a function of x, write V (J = v o + W * where v 0 ' 
projection of doublet velocity on the horizontal plane. Write 


% • ‘o 1 * v * % k 


(*.«) 


where t , m , n are the directional cosine of the unit vector n . Form the 
ooo ° 


vector quantity S by demanding S = V 1 x W 
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5 = r ( -ll - r Q n sin(fiT)- c/4n cos(nx)cos a Q )1 + (r Q n cos At 


cQ/4 sin(flr)cos a)j] * Wk 


(4.9) 


= W(U + r n sin(OT)+ c/4n cos(£Jt)cos a )j + W(r n cos(nx) 
0 0 0 


- cfl/4 sin(n-r)cos a)1 
It is now required that 


n Is 
u 


and 


♦ ♦ 

n 1 V n 
o o 


(4.m) 


n n = 1 


From (7), (R) and (9), there is obtained 


Wfil + r Q n sin(nT)+ c/4n cos(nx)cos a Q ) 


|/p ^ ^ 

I \r v + v 

r o o 


m = - 
o 


w(r Q n cos nt - cn/4 sin nx cos a Q 


j/u 2 * 


77 


P * v 


V 0 ' is the nagnitude of V Q ' 


(4.11) 


where 
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V ' 
o 


2 


(II f r o Q s1n(ftr)+ c/4a cos(nT)cos a Q ) * (r Q fl cos(n-r) 
- cfl/4 sinffiT)cos a Q ) ? (4.1?) 


Ry the sane procedure, n is ( rived by naklng the appropriate substitution , 

♦ ♦ ♦ ♦ 

c + -c, a Q ♦ a, t ♦ t, etc. (i.e., n = £i + m j + nk) to get 


i . w(j + rfl sin(nt)- c/4n cos(nt)cos a) 

I/— T4 
yw c v »• v 

m _ -W(rrt cos(ftt)+ cfl/d sin(nt)cos a) 

i/T-TT- TT 

l/H V ♦ " 


(4.13) 


|/w ? + v' 


where 

. 9 


(ll + rfl sin(nt) - c/4ft cos(nt)cos a)‘ + (rn cos(nt.) + cn/4 si n(iit)cosa)^ 

(4.14) 


since 


♦ ♦ ♦ ♦ 

r = ri + yj + zk. 


♦ ♦ ♦ ♦ 
P = xi+yj+zk 
0000 


(«.15) 


a vector D is defined as P 


♦ ♦ ■* 

P - P 0 ; then the magnitude |p| in eq. (1) can 


be expressed as 
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|fi| * { [t)( T-t ) + r cos(Ut)- r Q cos(nxW c/4(s1n(flt)cos a + s1n(nx)cos a Q )^ 

2 

+ [r sin(flt)- r o s1n(flx)-c/4(cos (flt)cos a + cos(flx)cos o^)] 


+ CW(t-T) -c/4(s1n a 0 + sin a )J } 
The final equation nay now be written as 


?il/2 




n.n 0 -3(n.n n ( v n n ) 

W 


" • T-l I q( r .t) [ " /' u " 1 dxdr 

4up ' J -« o* ' 1 :*i 3 1 o 


(4.17) 


where rip is a unit vector in the ?> direction, namely, n n = 


This integral equation was solved for the unknown q(r 0 ,T) by using a 
collocation process based on the vortex lattice assumption (ref. 6), The kernel 
is singular when P = n, and this was handled by use of the finite part 
technique (ref. 5). The remainder of the integrations was accomplished 
r-'meral ly. 

5. RESULTS ANP PISCUSSIPN 

The foregoing procedure was applied to a rotor ha vi .ig the dimension shown 
in fig. 1. The rotor was divided into 5 equal spanwise sections. The vortex 
lattice approximation was used which resulted in the placement of a vortex at 
the 1/4 chi rd positon on each section and a downward point was located at. the 
center of each span section and was placed at the 3/4 chord position as 
indicated on the figure. 
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The esults of calculations are given in fig. 2-4 for three cases. For all 
calculation the rotor speed was held constant at 30 rad/sec and the vertical 
velocity was 20 ft /sec. On figure 2 is shown the overall lift curve slope, 

Cl q , for 100 ft/sec forward velocity plotted against azimuthal angle. For 
the range of angles from approximately 90 degrees to 270 degrees, the lift curve 
slope is nearly Cjnstant. However, for the position of the angular position 
between 31 S degrees to 67.5° the lift curve slope is considerably ’ower. It 
could be surmised that with the rotor in a trailing position, that the effective 
aspect ratio of the rotor is lower than at other azimuthal positions and would 
accordingly have a lower left curve slope. 

On figure 3, the forward velocity has been increased to 200 ft/sec. Note 
that a large increase in C|_ q is found at 225 degrees. This one point is 
suspect, in as much as the stall region is being approached and the theory is 
based on linear approximations. Note that the region close to 0 deg. azimuth 
Ci_ a > is again low, where the minimum Cl q is 3.4 as compared to 3.9 for 
the 100 ft/sec case. 

On figure 4, the case f or a forward velocity of 300 ft/sec is shown. For 
this case, the rotor is definitely stalled in the azimuthal region 9 = 270° 
to 315°. 

fi. CONCLUDING REMARKS 

A lifting surface theory has been developed for a helicopter rotor in 
forward flight for incompressible flow. The method utilizes the concept of the 
linearized acceleration potential, and makes use of the vortex lattice procedure 
which has been successfully used for the non-rotating case. 

Results are shown for three forward velocities, ll = 100, 200 and 300 
ft/sec. The vertical velocity was held constant at W = 20 ft/sec as well as the 
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rotation speed, n = 3d rad/sec. For regions when the rotor blade was in a 

stalled condition the analytical results are not applicable and these regions 

are easily discernible in the plots. A significant finding in the results is 

the degradation in the lift curve slope, C L ^ t i n the neighborhood of 0 = 

fl deg. It should be pointed out that if two-dimensional coefficients were used 

in an analysis that this degradation in Cl would not be found. 

a 
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Fig. 2 Lift Coefficient Slope Versus Azimuthal Angle for 
U = 100 ft/sec, W = 20 ft/sec, Q = 30 rad/sec 
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Fig. 3 Lift Coefficient Slope versus Azimuthal Angle for 
U = 200 ft/sec, W = 20 ft/sec, fi = 30 rad/sec 
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Fig. 4 Lift Coefficient Slope Versus Azimuthal Angle for 
U = 300 ft/sec, W = 20 ft/sec, = 30 rad/sec 


